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This paper studies the hnearized gravitational field in the presence of boundaries. For 
this purpose, ^-function regularization is used to perform the mode-by-mode evaluation 
of BRST-invariant Faddeev-Popov amplitudes in the case of fiat Euclidean four-space 
bounded by a three-sphere. On choosing the de Bonder gauge-averaging term, the re- 
sulting (^(0) value is found to agree with the space-time covariant calculation of the same 
amplitudes, which relies on the recently corrected geometric formulas for the asymptotic 
heat kernel in the case of mixed boundary conditions. Two sets of mixed boundary condi- 
tions for Euclidean quantum gravity are then compared in detail. The analysis proves that 
one cannot restrict the path-integral measure to transverse-traceless perturbations. By 
contrast, gauge-invariant amplitudes are only obtained on considering from the beginning 
all perturbative modes of the gravitational field, jointly with ghost modes. 
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One-Loop Amplitudes in Euclidean Quantum Gravity 
I. INTRODUCTION 



In the Euclidean functional integral approach to quantum gravity [1], one deals with 
amplitudes written formally as path integrals over all Riemannian four-geometries match- 
ing the boundary data on (compact) Riemannian three-geometries ^Si, hi^ and ^E2, h2^- 

To take into account the gauge freedom of the theory, the path-integral measure also in- 
cludes suitable ghost fields, described geometrically by a one-form, hereafter denoted by 

(f = (fi^dx^^ (see appendix), subject to boundary conditions at (T>i,hij and ^S2,/i2^- Al- 
though a rigorous definition of the Feynman sum [2] over all Riemannian four-geometries 
with their topologies does not yet exist, the choice of boundary conditions still plays a key 
role to obtain a well-defined elliptic boundary-value problem, which may be applied to the 
semiclassical analysis of the quantum theory. 

In quantum cosmology, it was proposed in Refs. [3,4] that no boundary conditions 

should be imposed at the three-geometry (T,i,hi^, since this might shrink to a point in 
the case of the quantum state of the universe. One would then have to impose suitable 
boundary conditions only at ^^2, h2^, by describing the quantum state of the universe in 
terms of a Euclidean path integral over all compact Riemannian four-geometries match- 
ing the boundary data at ^E2,/i2^- Although this approach to quantum cosmology still 

involves a number of formal definitions, the semiclassical evaluation of the corresponding 
wave function may be put on solid grounds. The one-loop analysis is related to math- 
ematical and physical subjects such as cobordism theory (i.e. under which conditions a 
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compact manifold is the boundary of another compact manifold) , the geometry of compact 
Riemannian four-manifolds, the asymptotic heat kernel, the one-loop effective action and 
the use of mixed boundary conditions in quantum field theory (see below). 

In particular, over the last ten years many efforts have been produced to evaluate 
one-loop quantum amplitudes for gauge fields and the gravitational field in the presence 
of boundaries, either by using the space-time covariant Schwinger-DeWitt method [5] or 
the mode-by-mode analysis wich relies on ^-function regularization [6,7]. 

The main motivations were the need to understand the relation between different ap- 
proaches to quantum field theories in the presence of boundaries, and the quantization 
of closed cosmologies. Indeed, boundaries play an important role in the Feynman path- 
integral approach to quantum gravity [8] as we just said, in choosing BRST-covariant and 
gauge-invariant boundary conditions for quantum cosmology [9,10] and in studying differ- 
ent quantizazion and regularization techniques in field theory. In particular, for the latter 
problem, discrepancies were found in the semiclassical evaluation of quantum amplitudes 
by using space-time covariant methods, where the scaling factor of one-loop quantum am- 
plitudes coincides with the Schwinger-DeWitt A2 coefficient in the heat-kernel expansion, 
or noncovariant methods, where the same factor can be computed within the framework 
of ^-function regularization, and it is expressed through the ({0) value. 

If one reduces a field theory with first-class constraints to its physical degrees of 
freedom before quantization [11-14], one of the main problems is whether the resulting 
quantum theory is equivalent to the theories relying on the Faddeev-Popov gauge-averaging 
method or on the extended-phase-space Hamiltonian path integral of Batalin, Fradkin and 
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Vilkovisky, where one takes into account ghost and gauge modes [13-14]. We will see 
that, in a manifestly gauge-invariant formulation of such theories, there seem to be no 
unphysical modes, in that there are no subsets of the set of all perturbative modes whose 
effects cancel exactly the ones of ghost modes. This lack of cancellation turns out to be 
essential to achieve agreement between different techniques (see Sec. IV). 

In Ref. [11], the ({0) calculation was performed for gravitons by restricting the path- 
integral measure to transverse-traceless perturbations in the case of flat Euclidean four- 
space bounded by a three-sphere. In Refs. [15-17], this result was generalized to the 
part of the Riemannian de Sitter four-sphere bounded by a three-sphere. Both results did 
not coincide with those obtained by a space-time covariant method [9] . Hence the natural 
hypothesis arises that the possible non-cancellation of the contributions of gauge and ghost 
modes can be the cause of the discrepancy. In Refs. [18-21] such a suggestion was checked 
for the electromagnetic field on different manifolds and in different gauges. 

In Ref. [22] the asymptotic heat kernel for second-order elliptic operators was obtained 
in the case of pure and mixed boundary conditions in real Riemannian four- manifolds, and 
in Ref. [23] this analysis has been improved. In the light of these results, the conformal 
anomalies on Einstein spaces with boundaries have been re-calculated in Ref. [24]. 

In Ref. [25], we calculated the C(0) value for gravitons in the de Bonder gauge on the 
part of flat four-dimensional Euclidean space bounded by two concentric three-spheres, tak- 
ing into account the contribution of gauge modes and ghosts. The result was in agreement 
with the space-time covariant calculation. Hence we here investigate the linearized grav- 
itational field in the geometric framework of Ref. [11] (i.e. fiat Euclidean space bounded 
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by a three-sphere), and we compare the resulting ^(0) value with the space-time covariant 
calculation of the same Faddeev-Popov amplitudes, by using the recently corrected geo- 
metric formulas for the asymptotic heat kernel in the case of mixed boundary conditions 
[24]. For our purposes, we use the version of the ^-function tecnique [6] elaborated in 
Refs. [15-17]. Hence we write /„(M^) for the function occurring in the equation obeyed 
by the eigenvalues by virtue of boundary conditions, and d{n) for the degeneracy of the 
eigenvalues parametrized by the integer n. One then defines the function 

oo 

/(M2, s)^Y1 1^ fr^{M') . (1.1) 

n=no 

Such a function has an analytic continuation to the whole complex-s plane as a meromor- 
phic function, i.e. 

«/(m2, sY = + /^(M^) + 0(s) . (1.2) 

s 

Then the ({0) value is 

C(0) = /log + /pole(oo) - /pole(O) , (1.3) 

where Iiog is the coefficient of log M from I{M'^,s) as M — > oo, and Ipo\e{M'^) is the 
residue at s = 0. The uniform asymptotic expansions of basis functions as both their order 
and M tend to oo, yield /log and /poie(oo), while the limiting form of basis functions as 
M ^ and rz ^ oo yields /poie(O) [15-17,19-21,25]. 

In Sec. II we consider the mixed boundary conditions of Refs. [9,10] and compute 
the ^(0) value for the linearized gravitational field on the part of fiat four-dimensional 
Euclidean space bounded by a three-sphere, taking into account the contributions of gauge 
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modes and ghosts. Sec. Ill studies instead the mixed boundary conditions first proposed in 
Ref. [26] . Sec. IV compares the results of Sec. II with the ones deriving from a geometric 
analysis of the asymptotic heat kernel. Concluding remarks are presented in Sec. V, and 
relevant details are given in the appendix. 

II. LUCKOCK-MOSS-POLETTI 
BOUNDARY CONDITIONS 

In this section we evaluate ^(0) for the linearized gravitational field in the de Bonder 
gauge [25] 

^f'ih) = V^(V - ^^M-^^'V) > (2-1) 

where h is the perturbation of the background four-metric g, and V is the Levi-Civita 
connection compatible with g, i.e. Vg — 0. This gauge-averaging functional leads to 
the familiar form of the elliptic operator acting on metric perturbations. For example, in 
the case of a flat Euclidean background, this operator reduces to — [^ = — V^V^ [25]. 
Our background four-geometry is indeed flat Euclidean space bounded by a three-sphere, 
studied also in Ref. [11], and we take into account the contributions of gauge modes and 
ghosts. 

In Refs. [9,10] mixed boundary conditions have been introduced for quantum gravity 
and quantum cosmology. They are motivated by the need to obtain BRST-covariant 
boundary conditions which lead to gauge-invariant quantum amplitudes (see, however. 
Sec. V). This request leads to Dirichlet boundary conditions on the normal component 
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(po of the ghost one-form, to Robin conditions on the spatial components (pi of the ghost 
one-form and to mixed boundary conditions on the perturbative modes of the gravitational 
field. 

The Luckock-Moss-Poletti (also referred to as LMP) boundary conditions read {n 
being the normal to the boundary) [9,10,25] 



dM 



= 0, 



(2.2) 



h 



dM 



0, 



0, 



dM 



(2.3) 



(2.4) 



dM 



= 0, 



(2.5) 



dM 



= 0, 



(2.6) 



where K^i, and = 5^ — rif^'ny are the extrinsic-curvature tensor of the boundary and 
the tangential projection operator respectively. 

The mode-by-mode form of the LMP boundary conditions appears already in Ref. 
[25], and hence we here re-write them only for the spatial components (pi of the ghost, i.e. 



difi 2 

^ "f^ 

or r 



0, 



(2.7) 



dM 



to correct an unfortunate mistake in Ref. [25], where the numerical coefficient of ^ is 
incorrectly given (cf. Eqs. (3.3), (3. 12)- (3. 13) therein). That mistake did not affect the 
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C(0) value for ghosts, but it would affect the one-boundary calculations we are studying. 
On inserting the mode- by- mode expansion of the components fi{x, r) of the ghost one-form 
(see the appendix and Ref. [25]) into the boundary conditions (2.7) one obtains 

^--m„ = at dM, (2.8) 
or T 

^-lp^ = at dM. (2.9) 

Following the notation of the appendix and the technique used in Ref. [25] , one finds the 
following contributions to C(0): 



278 

C(0)transverse— traceless modes — 5 (^-lO) 



C(0)scalar modes — 18 gQ 2 ^ 6 180 ~ 45 ' C^'^^) 



C(0)partially decoupled modes — 2 17 — 19 , (2.12) 



_ ir _ 41 _ 2 _ 31 _ 1213 

C(U)vector modes ~ gQ ~ 3 ~ JgQ ~ QQ ' (2.1o) 



21 

C(0)decoupled vector mode — ^ > 

(2.14) 



_ . 119 899\ _ 314 

C(0)scalar ghost modes — ~2 | -|- j — , (2-15) 



_ / 19 209\__133 

S (Ujvector ghost modes — ~ { ^20 360 / ~ 90" ' v "/ 



5 

C(0)decoupled ghost mode — 77 • 

(2.17) 
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For example, in (2.11) we evaluate /log by using the uniform asymptotic expansions of 
modified Bessel fuctions. After eliminating fake roots of order {An — 1) as M — > one finds 
that the coefficient of InM is —{An +1). Hence one obtains 

oo 2 5 15 

^log = E y(-1 - = -2Ch(-3) + 18 + 2 = -60 + + 2 ■ ^^"^^^ 

n=3 

To evaluate /poie(oo)) the structure of the resulting n-dependent coefficient in the 
eigenvalue condition is 

{n^ - D 

(j^{n) = 12n ) , , (2.19) 
(n^ — A) 

2 

and it is easy to see that the term ^^\n.{aoo{n)) does not contribute as n — > oo (i.e. it has 
no coefficient of ^), and hence /poie(oo) = 0. 

Last, as M — > and n — >^ oo, the n-dependent coefficient which contributes to /poie(O) 

is 

^ ' ^ ' n3(n-2)(n + 2)2 ^ ' 

Thus, using Stirling's asymptotic expansion of the F-function [27], and after taking the 
coefficient of ^ as n — > oo in ^ln((To(n)), one finds 



/pole(0) = 9+^ + ^. (2.21) 



Finally, by virtue of (2.10)-(2.17) one gets the full ^(0) value for gravitons 



C(0) = -^. (2.22) 
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Hence the contributions of gauge and ghost modes do not cancel each other (cf. Ref. [11]) 
and our result, as we will show in Sec. IV, coincides with that obtained by using the 
covariant Schwinger-DeWitt technique on the part of flat Euclidean four-space bounded 
by a three-sphere. 

III. BARVINSKY BOUNDARY CONDITIONS 



The boundary conditions studied in Sec. II are not the only possible set of mixed 
boundary conditions for quantum gravity. As shown in Ref. [26], one can also set to 
zero at the boundary the gauge-averaging functional, the whole ghost one-form, and the 
perturbation of the induced three-metric. With the notation of Sec. II, after making a 
gauge transformation of the metric perturbation /i^i^ according to the law [11] 



(3.1) 



one finds in the de Bonder gauge (denoting by A the eigenvalues of the elliptic operator 



A 



(3.2) 



for any background with Ricci tensor R^xv In our flat Euclidean background, the Ricci 
tensor vanishes, and on making a 3-1-1 split of the de Bonder functional and of the 
ghost one-form (^^, the boundary conditions proposed in Ref. [26] read 



h 



dM 



h 



dM 



0, 



(3.3) 
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dM 



dM 



dM 



0, 



= 



dM 



^0 



dM 



= 0, 



= . 



dM 



(3.4) 
(3.5) 
(3.6) 
(3.7) 



Note that the vanishing of the whole ghost one-form at the boundary ensures the invariance 
of the boundary conditions (3.3) under the transformations (3.1) (see also Eq. (5.1)). At 
that stage, the only remaining set of boundary conditions on metric perturbations, whose 
invariance under (3.1) is again guaranteed by (3.6)-(3.7), is given by (3.4)-(3.5) by virtue of 
(3.2). In this respect, these boundary conditions are the natural generalization of magnetic 
boundary conditions for Euclidean Maxwell theory, where one sets to zero at the boundary 
the tangential components of the potential, the gauge-averaging functional, and hence the 
ghost zero-form [18-21]. The boundary conditions (3.3)-(3.7) were considered in Ref. [26] 
as part of the effort to understand the relation between the wave function of the universe 
and the effective action in quantum field theory. The loop expansion in quantum cosmology 
was then obtained after a thorough study of boundary conditions for the propagator [26] . 

In the light of (3.3), the boundary conditions (3.4)-(3.5) lead to mixed boundary 
conditions on the metric perturbations which take the form 



dhoo 6 

\- -h. 

OT T 



= 



dM 



dhoi 3 1 dhoo 
+ -hoi - 



dr T 



2 dx^ 



(3.8) 



(3.9) 



dM 
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To evaluate the scaling behaviour of the corresponding one-loop amplitudes, it is necessary 
to write down the mode-by-mode form of the boundary conditions (3.8)-(3.9), (3.3) and 
(3.6)-(3.7). They lead to (see Eqs. (A1)-(A13) of the appendix) 

dan 6 I den 2 ^ . f.^.^ 

— — -F-On ^— ^&n = at dM , (3.10) 

dr T dr 

^ + -&n-^^^^a„ = at dM, (3.11) 

(IT T Z 

^^ + -c„ = at dM , (3.12) 
dr T 

= at dM , (3.13) 

en = at aM, (3.14) 

/n = at aM , (3.15) 

/c„ = at aM , (3.16) 

= at dM , (3.17) 

mn = at dM , (3.18) 

p„ = at 5M . (3.19) 

On using the technique outlined in the introduction and applied also in Sec. II and in Ref. 
[25], the corresponding contributions to C(0) are found to be (see the appendix for details) 
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C(0)transverse— traceless modes — T7~ ' (3.20) 

4o 



_ 1 1 _ 764 

C(0) scalar modes ~ "'"^ ~ gQ ~ "'" ~ JgQ ~ "45" ' (3.21) 
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C(0)partially decoupled modes — 2 15 = 17 , (3.22) 



_ 11 2 31 _ 494 

C(0)vector modes ~ "'"^ ~ gQ ~ 3 ~ JgQ ~ ^5" ' (3.23) 



15 

C(0)decoupled vector mode — ^ j (3.24) 



si '-'/scalar ghost modes — -^1 ggQ j ~ 45 ' {o.Z,0) 



41 31 \ 77 

C(0)vector ghost modes = ~2 | — ~ / ~ 90 ' (3.26) 



5 

C(0)decoupled ghost mode — 2 ' (3.27) 

Hence the full C(O) for linearized gravity subject to the Barvinsky boundary conditions 
[26] is found to be 

241 

C(0) = -— . (3.28) 

Note that the result (3.28) differs from the one obtained in the previous section: ^0) = 
— As far as we can see, this property reflects two different physical situations described 
by two different sets of mixed boundary conditions (see also the comments in Sec. V). In 
the following section, the result (2.22) is checked by using geometric formulas for the 
asymptotic heat kernel. 
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IV. GEOMETRIC RESULTS ON THE HEAT KERNEL 



In this section we calculate by means of a geometric method the Schwinger-DeWitt 
coefficient A2 for the case of Luckock-Moss-Poletti boundary conditions. The result ob- 
tained coincides with the one found in Sec. II by a mode-by-mode analysis. The space-time 
covariant formulas for a second-order elliptic operator on a manifold with boundaries in 
the case of pure and mixed boundary conditions were obtained in Ref. [22] and corrected 
in Ref. [23]. We use them in the form presented in Ref. [24]. 

The Schwinger-DeWitt coefficient A2 for the elliptic operator 

where is a gauge derivative with curvature -F^i^, can be written as 

16%^ A2 = / b2dii + / C2dfi . (4.1) 

Jm JdM 

The volume coefficient &2 is well-known [5] , while surface terms depend upon the choice of 
boundary conditions. We use mixtures of Dirichlet and Robin boundary conditions, 

P_(/) = 0, + n'^V,)P+0 = , (4.2) 

where P± are projection operators [9,10,24]. 

The results can be expressed in terms of polynomials in the curvature tensor R^i^afS 
of the background four-manifold and in terms of the extrinsic-curvature tensor of the 
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boundary (hereafter R is the trace of the Ricci tensor, and K is the trace of K^^^)^ i.e. 



O 1 « 

q=^K^+ -jK^^K - ^KK^.K^^'^ + AKR 



(4.3) 



and 



(4.4) 



For Dirichlet boundary conditions [24] 



(4.5) 



while for Robin boundary conditions [24] 



1 



q + —p 



360" 45 



i(x-lii)A-+l„"V,(x-ifl) 



4 / 1 ^ 3 

-3(^-3^' 



(4.6) 



For mixed boundary conditions [23,24] 



C2 = Tr 



(4.7) 



where Latin indices run from 1 to 3, Greek indices run from to 3, and the stroke denotes 
(as in the rest of our paper) covariant differentiation tangentially with respect to the 
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three-dimensional Levi-Civita connection of the boundary [22] . The Luckock-Moss-Poletti 
boundary conditions for gravitons imply that [9] 

(i"+gravitons) = n^u, (2n"n^ - ^"'^j , (4.8) 



and 



V'gravitons — . (4-9) 



Similarly, for ghosts one has 



i'+ghosts) =5^^-nX. (4.10) 



and 

V'ghosts = -y • (4.11) 

Note that, on taking traces for gravitons in Eq. (4.7), one has to use the generalized 
Kronecker symbol 

Now inserting (4.8) and (4.9) into Eqs. (4.1), (4.6)-(4.7) one obtains for gravitons 



98 

A2 , =-— . (4.12) 

-^gravitons Q \ / 



Analogously, using the expressions (4.10)-(4.11) one finds for ghosts 



268 

^2,Host. = • (4.13) 
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These results are in full agreement with those in Sec. II. 

V. CONCLUSIONS 

Our paper has studied one-loop quantum gravity with boundaries in the limiting case 
of small three-geometries, which is relevant for quantum cosmology [11]. This leads to 
the analysis of flat Euclidean four-space bounded by a three-sphere, with mixed boundary 
conditions on metric perturbations. The results of our investigation are as follows. 

First, we have completed and improved the analysis of Luckock-Moss-Poletti mixed 
boundary conditions [9,10] considered in our earlier work [25]. These are motivated by 
the analysis of BRST transformations at the boundary, and are more relevant for su- 
persymmetric theories of gravitation [7,10]. The eight contributions to the full ^(0) for 
Faddeev-Popov amplitudes have been calculated in detail by means of a mode-by-mode 
analysis, and the resulting ^(0) value has been found to agree with the geometric theory 
of the asymptotic heat kernel. The latter relies on the results appearing in Refs. [22-24] . 

Second, we have performed a detailed analysis of the mixed boundary conditions for 
quantum gravity proposed by Barvinsky in Ref. [26]. Their main merit is the invariance 
under the transformations (3.1) on metric perturbations. The technique of ^-function 
regularization, jointly with a mode-by-mode analysis, has made it possible to obtain the 
full C(0) value as in Eq. (3.28). Indeed, this result differs from the one obtained in Sec. II. 
Thus, different boundary conditions lead to different semiclassical wavefunctions, with the 
exception of Euclidean Maxwell theory, where magnetic or electric boundary conditions 
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correspond to the same semiclassical amplitudes in the Lorentz gauge [20-21]. It should 
be emphasized that, in Sec. II, the invariance of the boundary conditions (2.3) under (3.1) 
leads to the Robin boundary conditions (2.6) on the ghost one-form. However, this implies 
that the boundary conditions (2.2) are not invariant under (3.1), since 

hij - hij = -t- Kij(^o . (5.1) 
In other words, the right-hand side of (5.1) does not vanish at dM if (2.6) holds, because it 



reduces to 



by virtue of (2.5). Indeed, in Ref. [9] the authors acknowledge that 



the boundary conditions (2.2)-(2.6) are not entirely BRST invariant, but they say, without 
explicit proof, that this does not affect the gauge invariance of the resulting quantum 
amplitudes. Hence we found it appropriate to consider also the boundary conditions of 
Sec. Ill, which are instead completely invariant under the transformations (3.1). On the 
other hand, the boundary conditions of Sec. II are motivated by self-adjointness theory 
[9] and are in agreement with the results on boundary conditions for one-forms, which 
should be mixed for gauge fields [23,28]. Hence the boundary conditions of Sec. Ill are 
less natural in this respect. 

Third, our analysis proves that, also in the one-boundary problem, one cannot restrict 
the path-integral measure to transverse-traceless perturbations. By contrast, ghost modes 
and the whole set of perturbative modes of the gravitational field are all necessary to 
obtain gauge-invariant (one-loop) amplitudes. 

It now remains to be seen what happens on considering arbitrary relativistic gauges 
and curved backgrounds. As one already knows from Euclidean Maxwell theory, it is not 
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possible to diagonalize the operator matrix of the problem for arbitrary gauge conditions 
[21]. In this respect, it appears interesting to characterize all relativistic gauge conditions 
which enable one to express the perturbative modes as linear combinations of Bessel func- 
tions [20,21,25]. Moreover, it is necessary to obtain geometric formulas for <^(0) in the case 
of Barvinsky boundary conditions studied in Sec. III. These boundary conditions are not 
naturally expressed in terms of projection operators, and hence the method of Refs. [22-24] 
cannot be applied. Nevertheless, work is in progress on this problem by the author of Refs. 
[23,28] and his collaborators. Last, but not least, the analysis of non-relativistic gauges 
for Euclidean quantum gravity within the Faddeev-Popov formalism, and the non-local 
nature of the one-loop effective action, deserve careful consideration. 

If one wants to relate our analysis to the Lorentzian theory, there is also the problem 
of interpreting the quantum state corresponding to the boundary conditions of Sees. II 
and III. Moreover, the impossibility to restrict the measure of the Euclidean path integral 
to transverse-traceless perturbations raises further interpretive issues for the Lorentzian 
theory, where such a reduction to physical degrees of freedom is instead quite natural. 
Thus, although the path- integral formulation of Euclidean quantum gravity [29] does not 
provide a mathematically consistent theory of the quantized gravitational field, the detailed 
calculations and the open problems presented in our paper seem to add evidence in favour 
of quantum cosmology having a deep influence on modern quantum field theory in four- 
dimensions [7,30]. 
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Following Ref. [25] , the metric perturbations are expanded on a family of three-spheres 
centred on the origin as 



APPENDIX 



oo 




(^1) 




oo 



(n2 - 1) 




(A2) 





20 



One-Loop Amplitudes in Euclidean Quantum Gravity 

OO 

+ E 



n=3 



/n(r) sl-^{x) + S^-^{x) 



{A3) 



where g(") (x) , S^""^ (x) , G^f (x) are scalar, transverse vector, and transverse-traceless ten- 
sor hyperspherical harmonics, respectively, on a unit three-sphere with metric Cy . The 
components ^po{x,t) and (fii{x,T) of the ghost one-form are expanded as in (Al)-(A2), 
providing one replaces an{T) , bnir) , Cn{T) by the modes ln{T),mn{T) and Pni'r) respec- 
tively [25]. Note that, strictly, our ghost one-form corresponds to a ghost ry^ and an 
anti-ghost ry^. They obey the same boundary conditions introduced for yj^ [10], and their 
contribution to ^(0) is obtained by using the multiplicative factor -2, as in Sees. II-IV. 

In the one-boundary problems studied in our paper, regular perturbative modes are 
only obtained by setting to zero the coefficients multiplying the modified Bessel functions 
Km since such functions are singular at the origin of flat Euclidean four-space. Hence one 
finds for scalar- type gravitational perturbations (cf. Ref. [25]) 



an{T) = ^ [71 In (Mr) + 73/n-2(MT) + 74/n+2(MT)] , 



{AA) 



hn{r) = l2ln{MT) + (n + l)73/n-2(MT) - (n - l)74/n+2(MT) 



dn{r) = T 



in-\-V\ (n — I') 

-f2ln{MT) + ) -(73/n_2(MT) -f- y—l^J^^^{MT) 



(n-2) 



(n + 2) 



(A5) 



(A6) 



en{r) = T 371 /n (Mr) - 2^2ln{MT) - 73/n-2 (Mr) - 74ln+2(MT) 



(A7) 
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The basis functions for vectorlike gravitational perturbations are 

Cn{r) = (Mr) + £2/n-i(Mr) , (^8) 



/n(r) = r 



^ -eiln+i (Mr) + -^—e2ln-i (Mr) 



(n + 2)^ "^^^ ^ (n-2) 
and the basis function for transverse-traceless symmetric tensor harmonics reads 



(A9) 



kn = a^TlniMr) . (AlO) 
Finally the basis functions for ghosts are 

Inir) = ^ [niln+iiMr) + Ks/n-ilMr)] , (All) 

mn(T) = -(n - l)KiIn+l{MT) + (n + l)AC2/n-l(Mr) , (A12) 

Pnir) = ^In{Mr) . (A13) 

The equations (A4)-(A13) hold both in Sec. II and in Sec. Ill, and we here focus on the 
latter application, since the former was treated in detail in Ref. [25]. 

On inserting the form of anir), bnir), dn{T),en{T) into the boundary conditions (3.10)- 
(3.11) and (3.13)-(3.14), one gets an eigenvalue condition for coupled scalar modes given 
by the vanishing of the determinant of a 4 x 4 matrix. As M — > 0, fake roots of order 4n 
are found to arise. In the calculation of /log, the factors \/M and compensate each 
other as M — > cxD, and hence /log is found to be 



°° 2 1 

V = E y(-4") = -2Ci?(-3) + 18 = -- + 18 . (A14) 



n=3 
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Moreover, as M — > oo and n — > oo, the n-dependent coefficient in the eigenvalue condition 
takes the form (cf. Eq. (2.19)) 



Poo(n) = 12n |^,_^j ■ (A15) 



Since B22^ [g an even function of n, the term ^hi(poo(^)) has no coefficient of ^ as 
77. — > oo, and hence /poie(oo) = 0. 

By contrast, as M — > and n — > oo, the n-dependent coefficient in the eigenvalue 
condition takes the form (cf. Eq. (2.20)) 



Thus, by virtue of Stirling's asymptotic expansion [27] 



In r(n) - ( n - ln(n) -n+ iln(27r) + 11-^1+ 0(n-^) , (A17) 



the /poie(O) value, which is the coefficient of ^ in the asymptotic expansion as n — > oo of 

2 

^ln(po(?^)), turns out to be 



/ . (0) = ^_1 + 1_1 = H. M18) 
poie^ ^ 180 6 3 6 180 ^ ^ 



The same technique yields the contributions (3.22)- (3.27) to the full C(0); bearing in 
mind that, for decoupled modes, the contribution to ^(0) is only given by the Iiog coefficient 
[15-17]. 
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